arXiv: 1507.06158v2 [hep-ph] 14 Oct 2015 


CNU-HEP-15-04 


Bounds on Higgs-Portal models from the LHC Higgs data 

Kingman Cheung 1,2,3 , P. Ko 4 , Jae Sik Lee 3,5 , and Po-Yan Tseng 1 
1 Department of Physics, National Tsing Hua University, 

101, Section 2 Kuang-Fu Road, Hsinchu, Taiwan 
2 Division of Quantum Phases and Devices, School of Physics, 

Konkuk University, Seoul If3-701, Republic of Korea 
3 Physics Division, National Center for Theoretical Sciences, 

101, Section 2 Kuang-Fu Road, Hsinchu, Taiwan 
4 School of Physics, KIAS, 85 Hoegiro Dongdaemun-gu, Seoul 02455, Korea 
5 Department of Physics, Chonnam National University, 

300 Yongbong-dong, Buk-gu, Gwangju, 500-757, Republic of Korea 

(Dated: October 15, 2015) 

Abstract 

In a number of Higgs-portal models, an SU( 2 ) isospin-singlet scalar boson generically appears 
at the electroweak scale and can mix with the Standard Model (SM) Higgs boson with a mixing 
angle a. This singlet scalar boson can have renormalizable couplings to a pair of dark matter 
particles, vector-like leptons or quarks, or new gauge bosons, thereby modifying the Higgs signal 
strengths in a nontrivial way. In this work, we perform global fits to such models using the most 
updated LHC Higgs-boson data and discuss the corresponding implications on Higgs-portal-type 
models. In particular we find that the current LHC Higgs-boson data slightly favors the SM over 
the Higgs-portal singlet-scalar models, which has to be further examined using the upcoming LHC 
Higgs-boson data. Finally, without non-SM particles contributing to the H 77 and Hgg vertices, 
the Higgs-portal models are constrained as follows: cos a > 0.86 and Ar tot < 1.24 MeV at 95 % 
confidence level (CL). 
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I. INTRODUCTION 


After the discovery of the 125 GeV boson at the LHC DEI, many analyses have been 
performed on the Higgs-boson data in order to identify the nature of the observe boson. So 
far, its properties are very close to those of the Standard Model (SM) Higgs boson: namely, 
(i) the spin, parity, and charge conjugation quantum numbers are equal to J PC = 0 ++ , 
which are in accordance with the SM Higgs boson, and (ii) its couplings to the SM particles 
are close to those of the SM Higgs boson at the end of the LHC Run I, which is indeed a 
remarkable achievement. The Higgs couplings to the SM particles are often parameterized 
in terms of the k’s defined as follows [3]: 

2 _ r i H -t a) __2 _ r to t(tf) + Ar tot ^ 

T^fTT • 1 -p > V D 

T(Lz —s- n)sM Tsm 

where i = W, Z, /, g, 7 , and Tsm denotes the SM total decay width while T tot (H) and Ar tot 
denote, respectively, the total decay width into the SM particles with modified couplings 
and an arbitrary non-SM contribution to the total decay width. The current best fits to the 
/q’s for i = W, Z, f from the ATLAS [1] and the CMS [5] collaborations are summarized in 
Table E 


TABLE I. The best fit values of k’s from the ATLAS Hj and CMS 0 at the end of the LHC Run 
I. The errors or the ranges are at 68% CL. 



KW 

Hz 

Ht 

K b 

K, r 

ATLAS 

0 .68±g;f° 

0-95t°;^ 

[- 0 . 80 , - 0 . 50 ] U [ 0 . 61 , 0 . 80 ] 

[-0.7, 0.7] 

[- 1 . 15 ,- 0 . 67 ] U [ 0 . 67 , 1 . 14 ] 

CMS 


bOS-oue 

o-si±8:ii 

0 74 + 0-33 
u - '+- 0.29 

0-84±8il 


New physics beyond the Standard Model (BSM) will be manifest itself if Ki 7 ^ 1 for some 
i in this approach. Very often it is assumed that the new physics effects are decoupled from 
the SM sector, thereby can be described by nonrenormalizable higher dimensional operators 
[ 6 j. This assumption encompasses a large class of BSMs, but still leaves out another large 
class of BSMs with an isospin-singlet scalar boson (of a mass around the electroweak (EW) 
scale) that could mix with the SM Higgs boson. This singlet scalar boson itself can couple 
to new particles such as a pair of dark matter (DM) particles, new vector-like quarks and/or 
leptons, new charged or neutral vector bosons, etc., just to name a few (see Ref. [7J for 
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more comprehensive discussion). Such a mixing between the singlet scalar boson and the 
SM Higgs boson does not decouple and cannot be captured by the usual higher dimensional 
operators, and therefore has to be treated in a separate manner. 

In Ref. [7j, a new parameterization was proposed which is suitable in the presence of a 
new singlet scalar boson that mixes with the SM Higgs boson. The singlet-mixed-in case 
deserves closer investigation, because many BSMs with good physics motivations come with 
an extra singlet scalar boson that can mix with the SM Higgs boson. This includes a large 
class of hidden-sector dark matter models such as Higgs-portal fermion or vector DM models, 
and DM models with local dark gauge symmetries, as well as nonsupersymmetric U(1)b-l 
model, vector-like fermions that could affect h —y gg, 77 , or models with the dilaton coupled 
to the trace of energy-momentum tensor. 

The Higgs-boson properties could be affected by the presence of new physics from different 
origins. The approach using /tj’s is simple and straightforward but in general it is difficult 
to further analyze the origin of new physics that had modified the s from the SM values. 
There are basically two different approaches to consider the new physics effects: one assumes 
either (i) the full SU(3)c x SU(2) L x U(l)y gauge symmetry or (ii) its unbroken subgroup 
SU(3)c x f/(l) em only in the effective Lagrangian for Higgs physics. Though either approach 
works as long as one is interested in any possible deviations of the Higgs couplings from the 
SM values, it would be more proper to impose the full gauge symmetry for investigations 
at the EW scale, because the energy and momentum transfer would be ~ 0(mz ) or higher. 
On the other hand, if we only impose the unbroken subgroup of the SM gauge group, 
the observed Higgs boson could be a mixture of the SM Higgs boson and other neutral 
scalar bosons that could mix with the SM Higgs boson after electroweak symmetry breaking 
(EWSB). Therefore, in order to isolate the effects of the mixing between the singlet scalar 
boson and the SM Higgs boson, we shall impose the full SM gauge symmetry when we 
construct the effective Lagrangian for the SM Higgs boson and the singlet scalar boson. 

In this paper, we perform the global fits to new physics scenarios with an extra singlet 
scalar boson mixed with the SM Higgs boson using the most recent Higgs data from LHC@7 
and 8 TeV. I 11 Sec. II, we set up the formalism used in this analysis, and compare it with 
the approach by the LHC Higgs Cross Section Working Group. In Sec. Ill, we give brief 
description of the models which are covered by our formalism. I 11 Sec. IV, we perform the 
numerical analysis with global fits to the LHC Higgs data, and present the best y 2 fit for 
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each model, and discuss the corresponding implications. Finally we summarize the results 
in Sec. V. 


II. FORMALISM 

In the following, we first describe the SM Higgs couplings to SM particles including 
fermions / and gauge bosons W, Z, 7 , g, and define a set of ratios bw,zj , 7 , s , which denote the 
size of the couplings relative to the corresponding SM one. Without loss of generality, we 
define a similar set of ratios cw,z,f,j,g for the singlet scalar boson couplings to the fermion / 
and gauge bosons W, Z, 7 , g relative to the corresponding one of the SM Higgs boson. After 
then we describe the mixing between the SM Higgs field and the singlet field via a mixing 
angle a. 


A. SM Higgs Couplings 

The couplings of the SM Higgs h to fermions are given by 

grrif 


hfl ~ ~ 


£h f 


2Mi 


w 


bfhff, 


( 2 ) 


f=u,d,l 

and its couplings to the the massive vector bosons by 


C hVV = gM w (b w W+W~» + b z 1 Z^zA h, 
V ^ 2 cos- Ow J 


(3) 


2 cos 2 9w 

where 6\y is the weak mixing angle. In the SM limit, we have bj = bw = bz = 1. 

While the SM Higgs coupling to two photons is defined through the amplitude for the 
decay process h —» 77 and it can be written as 


•A / 4'y'yh 


a M H <~, 7 / * * \ 

Anv bh (Clx ' e2j 


(4) 


where — 2ki(k 2 ■ ei )/Mj I , = e 2 — 2k 2 (k\ ■ e 2 )/M^ with ei ;2 being the wave 

vectors of the two photons and k\ ;2 being the momenta of the corresponding photons with 
( k\ + k 2 ) 2 = Mff. Including some additional loop contributions from non-SM particles and 
retaining only the dominant loop contributions from the third-generation fermions and W ± , 
the scalar form factor is given by 


Sl = 2 y NcQjbjF.f^-bwF^ + ASl = b,Si u , 

f=b,t,T 


(5) 
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where r x = Nc = 3 for quarks and Nq = 1 for taus, respectively. The additional 

contribution A S'l from non-SM particles is assumed to be real. Taking Mu = 125.5 GeV, we 
find S$ M = —6.64 + 0.0434 i. For the loop functions and the normalization of the amplitude, 
we refer to Ref. [ 8 j. 

The SM Higgs coupling to two gluons is given similarly as in h —> 77 . The amplitude for 
the decay process h —> gg can be written as 


■M-ggh — 


a,Mfj 5 ab 


47T V 


SI 


: i_l 


■ e 


2±J 


( 6 ) 


where a and b (a,b — 1 to 8 ) are indices of the eight SU(3) generators in the adjoint repre¬ 
sentation. Again, including some additional loop contributions from new non-SM particles, 
the scalar form factor is given by 


Si = E h U/M + AS? = b, S§ M . (7) 

f~b,t 

The additional contribution A Sj[ is assumed to be real. Taking M H = 125.5 GeV, we hnd 
5 | m = 0.651 + 0.0501 i. 

Finally, for the SM Higgs coupling to Z and 7, the amplitude for the decay process 
h —* Z(k\, 61 ) 7 (^ 2 , 62 ) can be written as 


■M-Z'yh — ~ 


a 

2 nv 


Sh 7 [h ■ k 2 el ■ e * 2 - h ■ 62 k 2 ■ el] 


( 8 ) 


where are the momenta of the Z boson and the photon (we note that 2ki-k 2 = Mjj— Mf), 
and €12 are their polarization vectors. The scalar form factor is given by 



= 2 y Q,N£m} 

f=t,b,r 


Il — 2 sin 2 9 W Q 
sin 9w cos 9w 


2 

1 b f F} 0) + Mf cot 9 W b w F w + A 


(9) 

( 10 ) 


The additional contribution ASjf 7 is assumed to be real. Taking Mu = 125.5 GeV, we hnd 
Ssm = — 11-0 + 0.0101 i. For the loop functions and the normalization of the amplitude, we 
refer to Ref. m- 


B. Couplings of the singlet scalar and the mixing 

The relative strength of the couplings of the singlet scalar boson s before mixing can be 
defined similarly in terms of a set of ratios <7 (i = /, W, Z, 7, g). Here the <7 parameterize 
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the couplings of s to the SM particles in a way similar to those of the SM Higgs boson h: 


£*ff - - E 


f=u,d,l 


9 m f 

2 Mi 


w 


CfSff , 


( 11 ) 


£ s vv — 9 M w (c w W+W M + c z 


9 - 

2 cos 2 0 W 


s, 


S] = 2 E N c 0 2 f c f F Bf (T f )-c w F 1 (T W ) + ASl c 7 5 s V 

f=b,t,r 

S g s = E + c 9 5f M , 

f=b,t 

s? = 2 EQ J 3 -2sm 2 ^Q 2 ^ F (o) + M | cot ^ = c Z7 ^ . 

^ T sin t/vn cos c/vc 

Since all the relative couplings Cj’s come from nonrenormalizable interactions between 
the singlet scalar s and the SM particles, except for the Higgs fields, one can simply assume 
that Cj’s are naturally suppressed by a heavy mass scale or a loop suppression factor: 


•‘ 0 ” + 


g 2, m 2 


or 


'‘ 0 ” + 


g 2 m 2 


(4t r) 2 M 2 ’ M 2 ’ 

On the other hand, the relative couplings V s of the SM Higgs boson with deviations com¬ 
ing from higher dimensional operators or additional particles running in the loop can be 
expressed as 


; ‘l” + 


g 2 m 2 


or “1” + 


g 2 m 2 


(4vr) 2 M 2 ’ M 2 ’ 

where M is the mass scale of a new particle that has been integrated out, and m is the 
external SM particles with m M, and g is a typical coupling of the SM particle and the 
heavy particle. Note that there would be extra loop suppression factors ( 1/(4 tt) 2 ) if the 

relevant operators are generated at one loop level. The sizes of s and Cj’s then set the 
stage for our numerical analysis. 

One further complication comes from the mixing between the SM Higgs held h and the 
singlet held s. The two mass eigenstates H X 2 are related to the interaction eigenstates by 
an 50(2) rotation: 


Hi = h cos a — s sin a ; H 2 — h sin a + s sin a , 


( 12 ) 


with cos a = c a and sin a = s a describing the mixing between the interaction eigenstates h 
(remnant of the SM Higgs doublet) and s (singlet). In this work, we are taking H\ = H for 
the 125 GeV boson discovered at the LHC and H 2 can be either heavier or lighter than H\. 
We are taking cos a > 0 without loss of generality. 
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Then, the relative couplings of the observed Higgs boson H to fermions /, gauge bosons 
W, Z, 7 , g are then given by 

bi cos a — c t sin a (i = f,W, Z, 7 , g) . (13) 

We observe that ki = ( biC a — C{S a ) 2 , see Eq. (JT|) . Note that the loop-induced Higgs decay 
with i = g, 7 can be modified by several different origins; (i) from scalar mixing denoted 
by a , (ii) from the singlet scalar couplings denoted by c 9i7 , especially when the singlet 
scalar couples to extra vector-like quarks and/or leptons or charged vector bosons, (iii) from 
modifications of the top and/or W boson couplings in the loop which are denoted by b t and 
bw , which arise from higher dimensional operators involving the SM Higgs doublet and the 
SM chiral fermions with the full SM gauge symmetry, (iv) from the couplings of the SM 
Higgs doublet to extra vector-like quarks and/or leptons or charged vector bosons, and (v) 
from some new physics effects that directly modify the couplings of the SM Higgs interaction 
eigenstate. 

The Ki parameterization is effective and simple but is highly degenerate, since different 
values of c a ,bi, Ci can lead to the same value of K t . It would be impossible to separate the 
true origin of new physics generating k / 1 in the k parameterization. 


C. Signal strength 


The theoretical signal strengths may be written as 


fl[V,V)~iHV) PIP), 


(14) 


where V = ggF. VBF. VII. ti ll denote the Higgs production mechanisms: gluon fusion 
(ggF), vector-boson fusion (VBF), and associated productions with a V = W/Z boson 
(VH) and top quarks (ttH) and V = 77 , ZZ , WW, bb, rf the decay channels. 

More explicitly, we are taking 

/2(ggF) = ( b g c a - c g s a ) 2 , /2(VBF) = fi(VH) = ( b v c a - c v s a ) 2 , p,(ttH) = ( b t c a - c t s a ) 2 

(15) 


with V = Z, W and 


fi<P) 


B(H ->• V) 
b(h sm —> V) 


(16) 
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with 


B(H ->• V) = 


r (H -> V) 


(i biC a - CiS Q ) 2 B(H SM £>) 


(17) 


Ftot(H) + Ar tot r tot (77)/r S M + Ar tot /rgM 

where % = 7 ,Z,W,b and r for V = 'f'y, ZZ,WW,bb and tt, respectively. Note that we 
introduce an arbitrary non-SM contribution Ar tot to the total decay width. Incidentally, 


T tot {H) becomes the SM total decay width Tsm when c a = 1 , bf — by — 1, A S^ ,9,z: = O^J 
and Ar to t = 0. For more details, we refer to Ref. ra¬ 


in. MODELS 

In a number of phenomenologically well motivated BSM models, there often appears a 
SM singlet scalar boson that can mix with the SM Higgs boson. Adding an extra singlet 
field to the SM is the simplest extension of the SM Higgs sector in terms of new degrees of 
freedom. A singlet scalar boson s does not affect the p parameter at tree level, and is not 
that strongly constrained by the electroweak precision tests (EWPT). It can also make the 
electroweak phase transition strongly first order [28], and enables us to consider electroweak 
baryogenesis if there are new sources of CP violation beyond the Kobayashi-Maskawa (KM) 
phase in the SM with three generations. Finally, if we imposed a new discrete Z 2 symmetry 
s —> —s, the singlet scalar s could make a good dark matter candidate [2fj]. This is the 
standard list for the rationales for considering a singlet scalar s. 

However, there are many more interesting scenarios where a singlet scalar appears in a 
natural way and plays many important roles. Let us list some examples, referring to Ref. [7] 
for more extensive discussion. 


A. Dark matter models with dark gauge symmetries and/or Higgs portals 

First of all, let us consider DM models where weak scale DM is stabilized by some 
spontaneously broken local dark gauge symmetries pi-12U] . This possibility is not that often 
considered seriously. However if we remind ourselves of the logic behind U( l) em gauge 
invariance, electric charge conservation, existence of massless photon and electron stability 
and non-observation of e —> 1 / 7 , one would realize immediately the same logic could be 
applied to the DM model building. One might think that this assumption may be too 

* We note 6 7 , Sj .z 7 = 1 when bf = by = 1 and A S'l’ 9,Zl = 0. 
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strong, since the lower bound on the DM lifetime is much weaker than that on the proton 
lifetime. This is in fact true, but this can be understood since proton is a composite particle, 
a bound state of 3 quarks with color gauge interaction, and baryon number violating operator 
in the SM is dim-6 or higher. Likewise longevity of DM might be due to some new strong 
interactions that make DM particle composite. Also, considering all the SM particles feel 
some gauge interactions, it would be natural to assume that the DM also may feel some 
gauge interactions (see Ref. |2JLJ for a recent review). 

In the case the dark matter particle is associated with some dark gauge symmetries, 
there would generically appear a dark Higgs boson after dark gauge symmetry breaking. 
The original dark Higgs <f> would be charged under some local dark gauge symmetry, but it 
is a singlet under the SM gauge group in the simplest setup. And after dark gauge symmetry 
breaking, there would be dark Higgs boson /i$, which would mix with the SM Higgs boson 
via the Higgs-portal interaction, 

_ d). 

A Higgs-portal coupling as small as A h$ ~ 10~ 6 can thermalize the hidden sector DM 
efficiently On the other hand, the effects of such a small coupling would be very difficult 
to observe at colliders. 

Also, the dark Higgs can stabilize the EW vacuum up to Planck scale, as well as it can 
modify the standard Higgs inflation scenario in such a way that a large tensor-to-scalar 
ratio r ~ (0.1) could be possible [22], which is independent of the precise values for the 
top quark and/or Higgs boson masses. Although the dark Higgs boson was introduced in 
order to break the dark gauge symmetry spontaneously, it has additional niceties in regard 
of cosmology in the context of the EW vacuum stability and the Higgs inflation assisted by 
Higgs-portal interaction. 

Even if we relax the assumption of the local dark gauge symmetry and consider more 
phenomenological Higgs-portal DM models, there will still appear a singlet scalar boson that 
can mix with the SM Higgs boson, if the Higgs-portal DM is a singlet Dirac fermion [23, 21J 
or a vector boson [HI [251 E6]. Also, it can play an important role in DM phenomenology. 
For example, one can easily accommodate the galactic center y-ray excess by DM pair 

t See, for example, Sec. Ill E and Fig. 5 (right panel) in Ref. [3D] for more details. 
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annihilation into a pair of dark Higgs bosons, followed by dark Higgs decays into the SM 
particles mmmm- 

Furthermore, there could be non-standard Higgs decays into a pair of lighter neutral 
scalar bosons (namely the dark Higgs boson) or a pair of dark gauge bosons, in addition to 
a pair of dark matter particles. In this case, the total decay width of the observed Higgs 
boson would receive additional contributions from the final states with dark matter, dark 
gauge bosons, or dark Higgs bosons, which are parameterized in terms of AT to t- Therefore, 
we will take the deviation Ar to t hi the total decay width of the 125 GeV Higgs boson as a 
free parameter when we perform global fits to the LHC data on the Higgs signal strengths. 

These classes of BSM models are phenomenologically very well motivated, and they have 
very significant impacts on the observed 125 GeV scalar boson. Therefore, it is very impor¬ 
tant to seek for a singlet scalar that can mix with the SM Higgs boson in all possible ways. 
The phenomenology associated with the observed Higgs boson measurements is straight¬ 
forward. The signal strengths of the 125GeV Higgs boson are suppressed from “1” in a 
universal manner, namely independent of production and decay channels. Moreover, the 
125 GeV Higgs couplings to SM fermions and weak gauge bosons are all suppressed by cos a 
relative to the SM values. One can also search for the heavier Higgs boson in this type of 
Higgs-portal models [32] . 

In summary, hidden-sector DM models are characterized by bi = 1 and q = 0 with a few 
simple implications: 

• Couplings to the SM fermions and gauge bosons are all suppressed by the factor cos a. 

• Decay Width: T(H —> V) = cos 2 aT SM (H — > V) and T tot (H) = cos 2 «T SM . Note that 
the total decay width of the Higgs boson, including the non-SM decay modes, is given 
by r tot (H) + AT tot • 

• Signal strengths: J1(V,V) ~ j2(V) /1(D) = cos 2 a .^ r ° t / rsM independently of the pro¬ 
duction mechanism V and the decay channel 'D. 

• Varying parameters: cos a and AT tot . 

In terms of two free parameters cos a and AT to t, we perform the y 2 minimization procedures 
on the LHC Higgs signal strength data in the next section. 
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B. Non-SUSY U(1)b~l extensions of the SM 


Another interesting example of Higgs-portal models is the nonsnpersymmetric U(1)b-l 
extension of the SM plus 3 RH neutrinos, which is anomaly free, so that no new colored or 
EW charged fermions are introduced: 

£ = £sm - V(H, <L) - Qa N j$N?Ni + Y N jHKN + Ac.) , (18) 

where the scalar potential V(H, $) is given by 

V(H, $) = + ^|#| 4 - A^|#| 2 |$| 2 + ^|<L| 4 . (19) 

Here the SM singlet scalar carries B — L charge “2”, and after B — L symmetry breaking 
from the nonzero VEV of <h, the resulting singlet scalar <p will mix with the SM Higgs field. 

If the B — L gauge boson Z' is light enough, the observed Higgs boson can decay into a 
pair of Z' bosons through the mixing between the SM Higgs boson and the U (1)e-n-charged 
singlet scalar 4>b-l, if this decay is kinematically allowed. The current bound on this model 
from the Drcll-Yan process is that )£ a few TeV, so that Qb-l ~ (a few) xlCT 3 or less, 
for this to happen. In this case, the Higgs phenomenology is described by two parameters, 
coscc and Ar tot , as in DM models with dark gauge symmetry and/or Higgs portals. 

C. Vector-like fermions for enhanced #(125) —y 77 and/or #(125) —> gg 

Right after the first candidate signature for the SM Higgs boson at the LHC, a number 
of groups considered new vector-like fermions (quarks or leptons) in order to explain the 
excessive signal strength in the H —>■ 77 decay channel. When considering vector-like 
fermions, one often has to introduce a singlet scalar field at renormalizable interaction level. 
Note that vector-like fermions can not directly couple to the SM Higgs doublet, and one 
has to introduce a singlet scalar coupled to them. Only in the presence of a singlet scalar, 
therefore, they can couple to the SM Higgs through the mixing between the SM Higgs boson 
and the new singlet scalar 0 In this case, the mixing between the SM Higgs boson and the 
new singlet scalar tends to reduce the signal strength for H 77 decay channel, although 
the loop contributions from the vector-like fermions would generate the singlet scalar decays 

t More detailed discussions of this class of models can be found in Sec. 3.4 in Ref. [7]. 
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into 77 and/or gg. It is essential to consider the mixing effects in the proper way (see, for 
example, Ref. [33]). 

In these types of models, the observed Higgs-boson couplings are given by 

gfjff = (bf cos a — Cf sin a) = cos a ; 

9hvv = (by cos a — cy sin a) = cos a for V = Z, W ; 

S]i 9 ^ = (S'h’ 9 ’ z ' r cos a - S]’ 9,Zl sin a) 

= cos a Rsm ’^ 7 + (AS'Z ’ 9 ’^ 7 cos a — A S]’ 9,Zl sin otj 
= cosaS^ + ASJi 9 ’ z \ (20) 

assuming bf = by = 1 and Cf = cy = 0 . 

In this case, the varying parameters are cos a , Aand/or A S 9 s , and possibly including 

Ar tot . 

D. Summary of the models 

Here we summarize the models in which a singlet scalar boson mixes with the SM Higgs 
boson: see Table [XT] for the relevant cf’s. More details including the corresponding La- 
grangian for each model can be found in Ref. [7J. 

Note that those classes of BSMs described in the subsections III.A and III.B are phe¬ 
nomenologically very well motivated by dark matter and neutrino physics as well as grand 
unification. Also, their impacts on the observed 125 GeV scalar boson as well as on the EW 
vacuum stability or Higgs inflation are straightforward: 

• The signal strengths of the 125 GeV Higgs boson are suppressed from “1” in a universal 
manner, namely independent of production and decay channels. 

• The 125 GeV Higgs couplings to the SM fermions and the weak gauge bosons are all 
suppressed by cos a relative to the SM values. 

• The additional singlet scalar boson can improve the stability of EW vacuum up to the 
Planck scale m- 

• The singlet scalar can improve the EW phase transition to be more strongly first order. 
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TABLE II. Nonvanishing cp ’s in various BSMs with an extra singlet scalar boson. For non-SUSY 
model, there would be nonzero c z > where Z' is the U(1)b-l gauge boson. Since the 
bound from Drell-Yan is very stringent, we will ignore H{ 125) — y Z Z' , although it is in principle 
possible if the gauge coupling is very small gs-L ~ a few xlCD 3 . Details can be found in Ref. [7J. 


Model 

Nonzero cp’s 

Pure Singlet Extension 

C/i 2 

Hidden Sector DM 

CpUi 2 

non-SUSY U(1) B -l 

c z ', c h 2 

Dilaton 

Cg 5 CW 5 CZ } Cry ? Cfo 2 

Vector-like Quarks 

Cgi C 7 j , Cfa2 

Vector-like Leptons 

c 1 ,cz 1 ,c h 2 

New Charged Vector bosons 

^75 Cfi2 

Extra charged scalar bosons 

Cg, C-y, CZ 7 , Chr 


• The mixing between the SM Higgs boson and the singlet scalar boson can modify the 
predictions for the tensor-to-scalar ratio within the Higgs inflation, and disconnecting 
the strong correlation of the inflationary observables from the top quark and the Higgs 
boson masses. 

Therefore it is very important to seek for a singlet scalar boson that can mix with the SM 
Higgs boson in all possible ways. 


IV. RESULTS 

We are going to perform the following fits: 

• SD fit - Singlet Dark Matter model and non-SUSY U(1)b-l case: Varying c a and 
AT tot , 

• SL fit - Singlet plus a vector-like Lepton: Varying s a , AT tot , A S], and AS], 
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TABLE III. The best-fitted values for SD, SL, and SQ fits. The SM chi-square per degree of 
freedom is Xg M /d.o.f.= 16.76/29, and p-value= 0.966. 


Fits 

X 2 

y 2 /dof p -value 


ATtot [MeV] 

AS h 

Best-fit values 

AS] AS] 

AS] 

Ca 

AS^ 


SD 

16.76 

0.621 

0.937 

0.000 

0.000 

— 

— 

- 

- 

1.000 

— 

- 

SL 

15.66 

0.626 

0.925 

0.129 

0.137 

-2.953 

-16.27 

- 

- 

0.992 

-0.835 

- 

SQ 

15.59 

0.678 

0.872 

0.036 

0.357 

0.875 

46.84 

1.315 

35.54 

0.999 

-0.832 

0.019 


• SQ fit - Singlet plus a vector-like Quark: Varying s a , Ar to t, AS 1 // AS], A S], and 
A SI. 

Note, instead of c a we vary s a in the SL and SQ fits because we have to specify c a and s a 
simultaneously in these fits. Otherwise, one may possibly explore the unphysical regions of 
ATtot < 0 and c a > 1 in the SD fit in order to study the parametric dependence. We neglect 
the Sh] couplings since we do not have any predictive power in the model-independent 
approach taken in this work. 

We use the most updated data summarized in Ref. [33] and the results of the fits are 
summarized in Table IH We find that the best-fit values of the SD fit are extremely close 
to the SM ones. For the SL and SQ fits, we observe that the best-fit values for AS]' 9 and 
AS]’ 9 are large while those for A S]j and ASfj are only about —0.8 and 0.02, respectively. 
In the remaining part of this Section, we discuss the details of each fit. 


A. SD 

In the SD fit, we scan the regions of parameters: c a C [0 : 2], Ar to t C [—4 : 8 MeV] 
including unphysical regions of Ar to t < 0 and c a > 1 to study the parametric dependence. 

In the left frame of Fig. [TJ we show the 68 % (Ay 2 = 2.30), 95 % (Ay 2 = 6.18), 99.7 
% (Ay 2 = 11.83) regions on the Ar tot -cosa: plane. When Ar to t > 0, we observe that the 
minima are developed along the yellow line in the black (Ay 2 < 0.01) region which is given 
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by the relation 


cos a = 


x l + i 1 + 4- 


Ar 


tot 


r 


SM 


1/2 


> 1 


( 21 ) 


In fact, the above relation can be obtained by requiring each signal strength to be the same 
as the SM one or jl(T,T>) = c^/(c 2 + Ar tot /r SM ) = 1. This implies the the best y 2 is 
obtained in the unphysical region of c a > 1. When c a < 1, we observe that the best y 2 is 
obtained again in the unphysical region of Ar to t < 0. If this is still the case in the future 
data, a large class of DM models (Higgs-portal fermion or vector DM, and DM models with 
local dark gauge symmetries) will be disfavored compared to the SM, except for the Higgs- 
portal scalar DM model without extra singlet scalar, for which the Higgs signal strength will 
be the same as the SM case. From the right frame of Fig. [TJ we see cos ct > 0.86 (0.81) and 
Ar tot < 1.24 (2) MeV at 95% (99.7%) CL. 


B. SL 

In the SL fit, we scan the regions of parameters: s a C [—1 : 1], Ar to t C [0 : 4 MeV], 
A S] C [-10 : 10], AS] C [-100 : 100], 

The CL regions are shown in Fig. [2j We observe that cos cc , sin a, Ar tot and AS# are 
well bounded as: 

cos a > 0.83 (0.76) at 95% (99.7%) CL ; 

| sin a | < 0.56(0.65) at 95% (99.7%) CL; 

Ar tot < 1.90 (3.00) MeV at 95% (99.7%) CL ; 

-2.95 (-3.96) < AS]j < 1.10 (2.02) at 95% (99.7%) CL . (22) 

In contrast, AS] and AS] are not bounded. From the relation AS# = AS] cosct — 
AS] since, in the limit since = 0, we see that AS] = AS# is bounded while AS] can 
take on any values. When | since| takes its largest value, AS] is most bounded: | AS'J | < 20 , 
see the lower-middle frame of Fig. [ 2 j As | AS'# | grows, a cancellation between the two terms 
AS] cosce and AS] since is needed to obtain the limited value of AS# together with non¬ 
vanishing AS] sin ce, explaining the wedges in the lower-left and lower-right frames. 

The best-fit values for AS] and AS] are —2.953 and —16.27, respectively, even though 
Ay 2 does not change much in most regions of the parameter space. Considering S] M = —6.64 
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and the best-fit value AS# = —0.835, let alone a certain level of cancellation, it would be 
very hard to achieve such large values for AS/ and AS/, unless the vector-like leptons are 
light, come in with a large multiplicity, and/or their Yukawa couplings to the singlet scalar 
s are strong. 

C. SQ 

In the SQ fit, we scan the regions of parameters: s a C [—1 : 1], Ar tot C [0 : 15 MeV], 
ASZ(ASJ) C [-10 : 10], ASJ(ASf) C [-100 : 100], 

The CL regions are shown in Fig. [3j We observe that cos cc , sin a, Ar to t, AS#, and 
AS#, are well bounded as: 

cos a > 0.70 (0.58) at 95% (99.7%) CL ; 

|sina| < 0.71(0.81) at 95% (99.7%) CL; 

Ar tot < 4.70 (10.40) MeV at 95% (99.7%) CL ; 

-2.94 (-3.96) < AS# < 1.14(2.05), at 95% (99.7%) CL; 

-0.13 (-0.18) < AS# < 0.35 (0.65) and 

-1.65 (-1.96) < A S 9 H < - 1.08 (-1.00) at 95% (99.7%) CL . 

One may make similar observations for AS/’ 9 and AS/’ 5 as in the SL case. The pa¬ 
rameters AS/’ 9 and AS/’ 9 are, in general, not bounded. When sin a = 0, AS/’ 9 = AS# 9 
and so they are bounded as AS# 9 . When | sina| takes its largest values, |AS/ ,9 | < 10. We 
also observe the wedges along sin a = 0 and AS/’ 9 = 0 due to the cancellation between 
AS/’ 9 cos a and AS/’ 9 since when |AS/’ 9 | > |AS# 9 |. 

The best-fit values for AS/(AS/) and AS/(AS 9 ) are 0.875(1.315) and 46.84(35.54), re¬ 
spectively, even though Ay 2 does not change much in most regions of the parameter space. 
Considering S/ M = —6.64 (Sf M = 0.65) and the best-fit value AS# = —0.832(AS# = 0.019), 
let alone a certain level of cancellation, it would be very difficult to achieve such large values 
for AS/’ 9 and AS/’ 9 , unless the vector-like quarks are light, come in with a large multiplicity, 
and/or their Yukawa couplings to the singlet scalar s are strong. 

Before closing this section, it is worth mentioning the experimental constraints on vector¬ 
like quarks and leptons. The vector-like quarks and leptons have been searched at the 
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Tevatron and at the LHC. The current best limits on vector-like quarks are from the ATLAS 
collaboration with 20.3 fb _1 luminosity at 8 TeV [53] , The limits range between 715 GeV 
and 950 GeV for up-type vector-like quarks and those for down-type ones between 575 GeV 
and 813 GeV. Considering these limits, we observe that one should have (9(100) vector-like 


quarks to accommodate the large best-fit values of AS]’ 9 shown in Table III assuming (9(1) 
Yukawa couplings of vector-like quarks to the singlet scalar s. 


V. DISCUSSION 

The Higgs-portal model involving a mixing between the SM Higgs field and an SU(2) 
singlet scalar boson is indeed the simplest extension to the SM Higgs sector, and gives rise 
to interesting phenomenology. In particular, this type of models can provide dark matter 
candidates, which exist in the hidden sector and interact with the SM sector through the 
mixing. Since it involves the mixing, so it will have non-negligible effects on the SM Higgs 
boson properties. In this work, we have used the most updated Higgs boson data from 
LHC@7 and 8 TeV to obtain very useful constraints on the models. In the simplest of this 
class of models - the singlet with a dark matter candidate (SD), the deviations from the 
SM Higgs couplings can be parameterized by the mixing cos a and the deviation in the total 
decay width AT tot . We found that the SD model does not provide a better fit than the SM, 
and thus we obtain the 95% CL on the parameters: 

cos a > 0.86 , AT tot < 1.24 MeV. (23) 

When more exotic particles are involved in the hidden sector, for example the vector-like 
leptons (SL) or vector-like quarks (SQ) in this work, the L /77 and Hgg vertices are modified 
non-trivially, and thus more parameters are involved. The constraints on cos cc and AT to t 
become somewhat less restrictive than the SD case (at 95%CL): 

SL : cosct > 0.83, AT to t ^ 1-9 MeV 
SQ : cos a > 0.70, AT to t ^ 4.7 MeV 

The allowed ranges for other parameters can be found in the previous section. 

We also offer the following comments on our findings: 

• The SM gives the best fit in terms of \ 2 /d.o.f. although the difference from other best 
fits (SD, SL, SQ) are not statistically significant yet. 
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• SD : In this case, the best y 2 occurs in the unphysical region: either c a > 1 or Ar < 0. 
If this is still the case in the future data, a large class of DM models (Higgs-portal 
fermion or vector DM, and DM models with local dark gauge symmetries) and non- 
SUSY U(1)b-l models will be strongly disfavored. However, the usual Higgs-portal 
scalar DM model with Z 2 symmetry without the extra singlet scalar may still be viable, 
since the Higgs signal strength in that model will be the same as the SM case. 

• SL : This case corresponds to the vector-like leptons in the loop for H —» 77 . We get 
a reasonably good fit. Nevertheless, we need a rather large value for AS] = —16.27, 
which might be possible only if the vector-like leptons are light, they come in with 
a large multiplicity, or the Yukawa couplings of the vector-like lepton to the singlet 
scalar s is strong. 

• SQ : This case corresponds to the vector-like quarks in the loop for H —> 77 and 
H —> gg. We get a reasonably good fit. However, we need a rather large value for 
AS] = 46.84, which might be possible only if the vector-like quarks are light, they 
come in with a large multiplicity, or their Yukawa coupling is very large. 

• SL and SQ : Though the best-fit values for AS]’ 9 and AS]’ 9 are large, those for 
A S]f = AS]' 9 cos a — AS]' 9 sin Q are only about — 0.8 and 0 . 02 , respectively. 

VI. ACKNOWLEDGMENTS 

The work of P.K. was supported in part by the Korea Neutrino Research Center (KNRC) 
which is established by the National Research Foundation of Korea (NRF) grant funded by 
the Korea government (MSIP) (No. 2009-0083526) and in part by the NRF Research Grant 
NRF-2015R1A2A1A05001869. The work of J.S.L. was supported by the National Research 
Foundation of Korea (NRF) grant (No. 2013R1A2A2A01015406). The work of K.C. was 
supported by the MoST of Taiwan under Grants No. NSC 102-2112-M-007-015-MY3. 


[1] G. Aad et al. [ATLAS Collaboration], “Observation of a new particle in the search for the 
Standard Model Higgs boson with the ATLAS detector at the LHC,” Phys. Lett. B 716, 1 
(2012) [arXiv:1207.7214 [hep-ex]]. 


18 



[2] S. Chatrchyan et al. [CMS Collaboration], “Observation of a new boson at a mass of 125 
GeV with the CMS experiment at the LHC,” Phys. Lett. B 716, 30 (2012) [arXiv: 1207.7235 
[hep-ex]]. 

[3] LHC Higgs Cross Section Working Group, A. David, A. Denner, M. Duehrssen, M. Grazzini, 
C. Grojean, G. Passarino and M. Schumacher et al., arXiv: 1209.0040 [hep-ph], 

[4] ATLAS Coll., “Updated coupling measurements of the Higgs boson with the ATLAS detector 
using up to 25 fb 1 of proton-proton collision data”, ATLAS-CONF-2014-009. 

[5] CMS Coll. (Vardan Khachatryan et al.), CMS-HIG-14-009, [arXiv: 1412.8662 [hep-ex]]. 

[6] A. Falkowski, arXiv: 1505.00046 [hep-ph]; J. Ellis, V. Sanz and T. You, JHEP 1407, 036 
(2014) [arXiv: 1404.3667 [hep-ph]]; E. Boos, V. Bunichev, M. Dubinin and Y. Kurihara, 
Phys. Rev. D 89, 035001 (2014) [arXiv: 1309.5410 [hep-ph]]; S. Banerjee, S. Mukhopadhyay 
and B. Mukhopadhyaya, Phys. Rev. D 89, no. 5, 053010 (2014) [arXiv: 1308.4860 [hep- 
ph]]; J. Elias-Miro, J. R. Espinosa, E. Masso and A. Pomarol, JHEP 1311, 066 (2013) 
[arXiv: 1308.1879 [hep-ph]]; B. Dumont, S. Fichet and G. von Gersdorff, JHEP 1307, 065 
(2013) [arXiv: 1304.3369 [hep-ph]]. 

[7] S. Choi, S. Jung and P. Ko, “Implications of LHC data on 125 GeV Higgs-like boson for the 
Standard Model and its various extensions,” JHEP 1310 (2013) 225 [arXiv: 1307.3948]. 

[8] J. S. Lee, A. Pilaftsis, M. S. Carena, S. Y. Choi, M. Drees, J. R. Ellis and C. E. M. Wagner, 
“CPsuperH: A Computational tool for Higgs phenomenology in the minimal supersymmetric 
standard model with explicit CP violation,” Comput. Phys. Commun. 156 (2004) 283 [hep- 
ph/0307377], 

[9] T. Hur, D. W. Jung, P. Ko and J. Y. Lee, “Electroweak symmetry breaking and cold dark mat¬ 
ter from strongly interacting hidden sector,” Phys. Lett. B 696, 262 (2011) [arXiv:0709.1218 
[hep-ph]]. 

[10] T. Hur and P. Ko, “Scale invariant extension of the standard model with strongly interacting 
hidden sector,” Phys. Rev. Lett. 106, 141802 (2011) [arXiv: 1103.2571 [hep-ph]]. 

[11] S. Baek, P. Ko, W. I. Park and E. Senaha, “Higgs Portal Vector Dark Matter : Revisited,” 
JHEP 1305, 036 (2013) [arXiv:1212.2131 [hep-ph]]. 

[12] S. Baek, P. Ko and W. I. Park, “Singlet Portal Extensions of the Standard Seesaw Models to a 
Dark Sector with Local Dark Symmetry,” JHEP 1307, 013 (2013) [arXiv: 1303.4280 [hep-ph]]. 


19 



[13] S. Baek, P. Ko and W. I. Park, “Hidden sector monopole, vector dark matter and dark 
radiation with Higgs portal,” JCAP 1410, no. 10, 067 (2014) [arXiv: 1311.1035 [hep-ph]]. 

[14] P. Ko and Y. Tang, “Self-interacting scalar dark matter with local Z 3 symmetry,” JCAP 
1405, 047 (2014) [arXiv: 1402.6449 [hep-ph], arXiv: 1402.6449]. 

[15] P. Ko and Y. Tang, VAMDM: A Model for Sterile Neutrino and Dark Matter Reconciles 
Cosmological and Neutrino Oscillation Data after BICEP2,” Phys. Lett. B 739, 62 (2014) 
[arXiv: 1404.0236 [hep-ph]]. 

[16] S. Baek, P. Ko and W. I. Park, “Local Z -2 scalar dark matter model confronting galactic 
GeV-scale 7 -ray,” Phys. Lett. B 747, 255 (2015) [arXiv: 1407.6588 [hep-ph]]. 

[17] P. Ko and Y. Tang, “AMS02 positron excess from decaying fermion DM with local dark gauge 
symmetry,” Phys. Lett. B 741, 284 (2015) [arXiv:1410.7657 [hep-ph]]. 

[18] P. Ko, W. I. Park and Y. Tang, “Higgs portal vector dark matter for GeV scale 7 -ray excess 
from galactic center,” JCAP 1409, 013 (2014) [arXiv: 1404.5257 [hep-ph]]. 

[19] P. Ko and Y. Tang, “Galactic center 7 -ray excess in hidden sector DM models with dark gauge 
symmetries: local Z 3 symmetry as an example,” JCAP 1501, 023 (2015) [arXiv: 1407.5492 
[hep-ph]]. 

[20] P. Ko and Y. Tang, “Dark Higgs Channel for FERMI GeV 7 -ray Excess,” arXiv: 1504.03908 
[hep-ph]. 

[21] P. Ko, “Dark matter, dark radiation and Higgs phenomenology in the hidden sector DM 
models,” arXiv: 1503.05412 [hep-ph], 

[22] P. Ko and W. I. Park, “Higgs-portal assisted Higgs inflation with a large tensor-to-scalar 
ratio,” arXiv: 1405.1635 [hep-ph], 

[23] S. Baek, P. Ko and W. I. Park, “Search for the Higgs portal to a singlet fermionic dark matter 
at the LHC,” JHEP 1202, 047 (2012) [arXiv: 1112.1847 [hep-ph]]. 

[24] S. Baek, P. Ko, W. I. Park and E. Senaha, “Vacuum structure and stability of a singlet fermion 
dark matter model with a singlet scalar messenger,” JHEP 1211, 116 (2012) [arXiv: 1209.4163 
[hep-ph]]. 

[25] Y. Farzan and A. R. Akbarieh, JCAP 1210, 026 (2012) [arXiv: 1207.4272 [hep-ph]]. 

[26] M. Duch, B. Grzadkowski and M. McGarrie, arXiv: 1506.08805 [hep-ph], 

[27] K. Cheung, J. S. Lee and P. Y. Tseng, “Higgs Precision (Higgcision) Era begins,” JHEP 1305 
(2013) 134 [arXiv:1302.3794 [hep-ph]]. 


20 



[28] D. A. Kirzhnits and A. D. Linde, Annals Phys. 101, 195 (1976); D. E. Morrissey and 
M. J. Ramsey-Musolf, New J. Phys. 14, 125003 (2012) [arXiv: 1206.2942 [hep-ph]]. 

[29] V. Silveira and A. Zee, Phys. Lett. B 161, 136 (1985); C. P. Burgess, M. Pospelov and T. ter 
Veldhuis, Nucl. Phys. B 619, 709 (2001) [hep-ph/0011335]; J. M. Cline, K. Kainulainen, 
P. Scott and C. Weniger, Phys. Rev. D 88 , 055025 (2013) [Phys. Rev. D 92, no. 3, 039906 
(2015)] [arXiv: 1306.4710 [hep-ph]]. 

[30] S. Baek, P. Ko, W. I. Park and Y. Tang, JCAP 1406, 046 (2014) [arXiv:1402.2115 [hep-ph]]. 

[31] C. Boehm, M. J. Dolan and C. McCabe, “A weighty interpretation of the Galactic Centre 
excess,” Phys. Rev. D 90 (2014) 2, 023531 [arXiv: 1404.4977 [hep-ph]]. 

[32] S. Banerjee, M. Mitra and M. Spannowsky, arXiv: 1506.06415 [hep-ph]; D. Buttazzo, F. Sala 
and A. Tesi, arXiv: 1505.05488 [hep-ph]. 

[33] T. Abe, N. Chen and H. J. He, “LHC Higgs Signatures from Extended Electroweak Gauge 
Symmetry,” JHEP 1301, 082 (2013) [arXiv:1207.4103 [hep-ph]]. 

[34] K. Cheung, J. S. Lee and P. Y. Tseng, “Higgs precision analysis updates 2014,” Phys. Rev. 
D 90 (2014) 095009 [arXiv: 1407.8236 [hep-ph]]. 

[35] G. Aad et al. [ATLAS Collaboration], “Search for production of vector-like quark pairs and 
of four top quarks in the lepton-plus-jets final state in pp collisions at y/s = 8 TeV with the 
ATLAS detector,” JHEP 1508 (2015) 105 [arXiv: 1505.04306 [hep-ex]]. 


21 




FIG. 1. The 68 % (Ay 2 = 2.30), 95 % (Ax 2 = 6.18), 99.7 % (Ax 2 = 11.83) confidence- 
level (CL) regions for the SD fit on the AFtot-cosa plane. The horizontal line in the left frame 
shows the physical limit cos a = 1. In the right frame, we show the CL regions after applying 
ATtat > 0 and cos a < 1. The best-fit points are along the yellow line passing through the point 
(AT, cos ct) = (0,1) in the left frame. In the black regions, we have Ay 2 < 0.01. 







FIG. 2. The CL regions for the SL fit. The description of the CL regions is the same as in Fig. [lj 
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